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Abstract 

We show that the effective running coupling constant, a e ff, and the effective regulator, Ktn 2 D , 
which we used recently to calculate the energy loss, and the elliptic flow, \>i, of heavy quarks 
in an expanding quark gluon plasma plasma (QGP) fit] are compatible with lattice results and 
with recently advanced analytical pQCD calculation. 



1. Introduction 

The analysis of the spectra of heavy mesons, observed by single non photonic electrons in 
heavy ion collisions at 200 AGeV QQ], have revealed two surprising results: a) Small values 
of Raa = dcr AA /dpj/(< N c > d<T pp /dpj), where < N c > is the average number of initial binary 
collisions b) A substantial elliptic flow x>2 =< cos 2(<p - <p re action) >■ Initially heavy quarks have 
no elliptic flow and without any interaction we expect Raa to be one (assuming that the parton 
distribution function is similar in pp and in heavy ions). The finite value of x>2 and the low 
Raa value point therefore to a substantial interaction with light partons during the expansion of 
the system. If in heavy ion collisions a QGP is formed the natural starting point to calculate 
the interaction with the constituents of the plasma, the light quarks and gluons, is perturbative 
QCD (pQCD). Comparison with data on Raa and on V2 shows that the interaction has to be 
much stronger than obtained in pQCD calculations using a standard fixed value of the coupling 
constant 0.3 < a s < .5 and a infrared regulator yt 2 in the order of the thermal squared Debye 
mass, m 2 D = (1 + ^■)4-na s T 2 , where Nf is the number of flavors. 

The weak point of several models aiming at reproducing these Raa and V2 observables is the 
assumptions on a s and on /i, both being not well determined and more or less taken as input 
parameters. Recently we advanced an approach in which we related these parameters to physical 
observables (U: 

• We employ a running nonperturbative effective coupling constant a e g(Q 2 ) whose value 
remains finite if the Mandelstam variable t approaches by adopting the parametrization of 
Ref. [J] for the time-like sector and by truncating the 1-loop renormalized coupling constant in 
the space-like sector in order to satisfy the so-called "universality constrain" postulated in 0], 
which leads to 

4 ( i-^'atn(L + /7r) for Q 2 > 

« -» «eff(e 2 ) = — a S a, for - |Q 2 |sat < Q 2 < (1) 

Pq 1 LZ 1 for Q 2 < -\Q 2 U 
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with A, = 11 - | N f , N f = 3, ff sat = 1.12, ie 2 | sat = 0.14GeV 2 andL ± = \n(+Q 2 /A 2 ) . 

• We extend the usual Braaten-Yuan scheme, applied by Braaten and Thoma |7] (BT) to 
fixed-a, QCD, to the case of a running coupling constant, introducing a semi-hard propagator 
with m 2 eff (T, Q 2 ) = (l + Ana eS {Q 2 )T 2 , in the \Q 2 \ > \t* | sector in order to guaran- 



tee a maximal independence of the energy loss *P w.r.t. this unphysical t* scale. 

• We then define an effective one massive-gluon exchange model (OGE) by replacing the 

standard polarization in the t-channel [6] /j 2 = where mo is the Debye mass, by an effective 
one, i.e. p? = kiii 2 d , where ih 2 D (T) — (\ + ^}4-na s (-m£,{T) 2 )T 2 is a self-consistent Debye 
mass, k is determined by calibrating ^| to the one calculated in our (extended) hard thermal loop 
calculation. The details of this procedure is found in the appendix of ref. flJJ] and the details of 
the OGE models labeled "model E" in the text of the same reference. 

With these improvements we were able to reproduce the centrality dependence of the exper- 
imental Raa up to a factor of 2-3 as well as the experimental value of V2 with collisional energy 
loss only. For details we refer to QJJ] . 

It is the purpose of this contribution to show that our choice of a e s and p are compatible with 
lattice data. Kaczmarek and Zantow [8] have studied the potential and coupling constant on the 
lattice in 2-flavor QCD by investigating the free energy between two heavy quarks. Assuming 



that at small distances the singlet potential has the form V{r) 



4 cc. T 



-f they define a qq {r, T) 



3 r 2£^i£) j n or( j er to compare our approach with these lattice data we start out from the t-channel 
matrix element, which consists of two parts: a hard thermal loop (HTL) approach for {\Q 2 \ < \t*\) 
and a Born "semi hard" amplitude for (\Q 2 \ > |f*|) 



V Pe (T;4) 



Ana eSS {-q ) 



&(q 2 < |f* |) 



Q(q 2 > |f* |) 



q 2 + fI L (w = 0, q) ql + A(p Q )m 2 D JT, -q 2 ) 



(2) 



with q = \\q\\ and Hl(o> — 0, q) = m 2 D eff (7\ -q 2 ) and |f*| chosen as T 2 . Notice that the parameter A, 
and thus our regulator, is found to depend on the heavy quark momentum pq in our prescription, 
as displayed on figure [I] (left), but neither on T nor on the heavy-quark mass M when displayed 
as a function of pg/Mj- For pq/M — > 0, A approaches 1 and our regulator comes into the range 
of the screening mass derived for a static potential in lattice QCD [8]. 

The potential in coordinate space is then obtained by Fourier transformation. The radial 
force - with is the essential quantity for 4| - we obtain by derivation. It is compared in fig[U 
(middle and right) with results from lattice calculations. We present our results for three values 
of A: A = (standard BT) and 1 (= pq = 0) present the boundaries where our approach is 
meaningful, A = 0. 1 1 (red full line) is the value we obtain by the above procedure for pq — > oo. 
The lattice results for the free energy F and for the energy U are the limits of the vertical lines. 
Whereas at very small distances as compared to the screening length the potential corresponds to 
the energy with increasing distance the polarization of the medium becomes important and the 
discussion whether the free energy or the energy presents the Schrodinger equivalent potential 
is not settled yet. For Pq = — strictly speaking, the only case where the lattice data are 
relevant for comparison - our procedure leads to A = 1 and a force close to the one deduced 
for the assumption V = F. For increasing pq, A decreases correlatively, and our force is in 
between the two lattice choices V — F and V = U up to r w 1 fm. Beyond r — 1 fm, and for 



Standard BT corresponds to A = 0, with large numerical dependences of ^ on t* . 
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Figure 1 : (Color online) Left: the A parameter (eqO which guarantees the maximal insensitivity of the momentum loss 
per unit length on the unphy sical scale f . Middle and right: Comparison of dV/dr for A = and 1 (full curves bordering 
the gray bands) and for A = 0.11 (thick full red line) and for model E (full dashed black line) of ref. [1] with lattice 
results (vertical lines whose limits correspond to the V = V — top - and V = F — bottom - prescriptions). 

T « 2T C we observe a rather long tail for small A's and thus relatively moderate values of the 
associate regulator. The physical relevance of this tail can only be addressed once lattice results 
for pq » M are available. An alternative way of comparing our model to lattice data is to resort 
to the effective T-dependent coupling constant a qq defined in Q8[]. It is presented in figH] 




Figure 2: (Color online) The running coupling constant a qq in our approach as compared to lattice results for the same 
approaches and conventins as in fig[T](middle and right). 



In the interesting regime around 0.5-1 fm, corresponding to 0.2 to 0.4 GeV momentum trans- 
fer, the coupling constant which we employ is around T c at the lower limit of that predicted by 
lattice QCD. Only around \.5T C we start to overpredict a s for r > 0.8 fm, but this region is of 
moderate (but finite) importance for the energy loss of the heavy quark. In figs. Q]and|2] the 
dashed lines correspond to the OGE model E introduced in ref. [1]. Let us remind that model 
E was obtained by fixing, for simplicity reasons, A = 0.1 1 for all values of pq and deducing an 
effective gluon square mass of Kin 2 D with k « 0.2. 

Our approach can also be compared with a recent analytical approach to derive a running 
coupling constant by Peigne and Peshier J9(]. In the limit of large plasma temperatures and 
therefore of a weak coupling constant as they obtain: 



dx T 



AnT 2 2 
= — — a s {m D )a s {ET) 

2 ET 
+ 9 Xn W 2+C(Nf) + ° 



l + N -A^ ET 



6 / " m 2 D 
ET^ 

a s (m 2 D ) In — T 



(3) 



with c(Nf) « 0.146./V/ + 0.05 and iho being the self-consistent Debye mass. Figj3]compares our 

3 




LOO 




50 



'-=4Ge\ 



20 
10 
5 



2 



5 10 



50 100 



5001000 



PQ[GeV] 



Figure 3: (Color online) Comparison of our energy loss with that of ref. |9j] which uses a running coupling constant 
calculated in pQCD. For the symbols we refer to the text. 



whereas our results are the light intervals (bounded by A — 0.1 1 and A - I). The full line is our 
result for A = 0.1 1. The dashed line represent the analytical result for a fixed coupling constant 
10. Thus in the range of validity of the approach of ref. jjj, at very high temperatures and/or 
Pq, our approach reproduces well the most recent pQCD calculations. 

In summary, we have shown that our choice of coupling constant and regulators, motivated by 
HTL-type calculations, agrees satisfactorily with lattice and analytical pQCD calculations. The 
associated drag force deduced in our approach is an increasing function of the temperature. In 
particular, it shows no peak around the transition temperature T c , contrarily to approaches where 
the qQ interaction is strictly taken as the internal energy J/EI. Phenomenological consequences 
of these patterns should be addressed systematically in the future. 
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